Trajectory-constrained optimal local time-continuous waveform controls for state 

transitions in iV-level quantum systems 
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Based on a parametrization of pure quantum states we explicitly construct a sequence of (at most) 
4N — 5 local time-continuous waveforms controls to achieve the specified state transition for iV-level 
quantum system when sufficient controls of the Hamiltonian are available. The control magnitudes 
are further optimized in terms of the time-energy performance which is the generalization of the 
time performance, and then the trajectory-constrained optimal local time-continuous waveforms 
controls including both local sine- waveforms and n-order-polynomial-function-waveform controls are 
obtained in terms of time-energy performance. It is demonstrated that constrained optimal local 
nst-order-polynomial-function-waveform controls approach constrained optimal bang-bang controls 
when n — > oo. 
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I. INTRODUCTION 



Control of quantum systems has been recognized as an 
important issue for some time with early beginnings of 
the application of control theory in the quantum domain 
dating back to the 1980s [3-l|. Generally, the main goal 
of control theory is to find controls leading the objects to 
a desired situation. There are usually two ways of speci- 
fying "a desired prescribed situation" : the controllability 
viewpoint and the optimization viewpoint Q. The con- 
trollability of quantum systems has been investigated by 
many researchers (see Ref. 0, ID, Q). Optimal control 
theory has also been successfully applied to the design of 
open-loop coherent control strategies in physical chem- 
istry 0, Q . Recently, time-optimal control problems for 
spin systems have been solved to achieve specified control 
objectives in minimum time [9Tjll|. In general, optimal 
control problems can only be solved by using numerical 
optimization techniques. 

When sufficient controls of the Hamiltonian are avail- 
able, the third kind of mixed approach may be adopted: 
one can first construct some simple controls along a 
chosen trajectory to achieve a desired state transition, 
and then make full use of degree of freedom to opti- 
mize some kind of performances including minimum time 
performance. In this way, trajectory-constrained sim- 
ple optimal controls are obtained. Present technology 
justifies this third approach for some super-conducting 
systems |12|- l l4 | and nuclear magnetic resonance sys- 
tems(NMR) [15|], for example. It has been demonstrated 
that simple waveforms such as local square wave func- 
tion (Bang Bang control) [21] can be constructed to 
achieve the desired state transition. Quite recently, Bang 



Bang control has been successfully applied in some phys- 
ical systems [lH- 20]. In practice the control functions 
need not be restricted to bang-bang controls but other 
time-continuous function waveforms such as triangle- 
waveforms and sine- waveform may be used to manipulate 
quantum states of N— level quantum systems. In this pa- 
per, the third kind of mixed viewpoint is adopted and it is 
demonstrated that local time-continuous function wave- 
forms including both local n-order-polynomial waveforms 
and sine waveforms can be constructed to achieve the de- 
sired state transitions. 



II. PREREQUISITE 

The state of an iV-level quantum mechanical system 
is presented by a vector in .ZV-dimensional Hilbert space 
H. In quantum mechanics, the state ip is denoted as 
and called ket. To any \ip) is associated a linear operator 
(ip\ : H — > C, which is called bra. Given a ket \tp) and a 
bra (</>], we define a linear operator |^)(^>| : H — s- H . The 
state of an iV-level quantum system can be expressed as 
= ^2n=a c n\ n ) with regard to a chosen basis |n) (n = 
0, 1,2,..., AT — l)in the Hilbert space. The coefficients c n 
are a set of complex numbers satisfying the normalization 
constraint J2n=o l c ™l 2 = By ignoring the global phase, 



the coefficients can be expressed as 
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with < 0i, ... , 9n-i < 7T and < <f>i, . . . , 4>n-i < 27r. 
Thus any pure state \ip) of an N— level quantum system 
can be represented by the 2(N — 1) generalized geometric 
parameters (0; $) with = (6%, . . . , 6n-i) T and $ = 
(</>i,...,</>jv_i) t . 
For N > 2 we define the operators (JV x iV matrices) 

x Ntk = \k){k+l\ + \k+l){k\ (2a) 

tfjv,fc = *[|*: + l)<A| -!*)<*: + 1|] (2b) 

z N ,k = lN-2\k+l){k + l\ (2c) 

7iv,fe = |fc)<fc| + |fe + l)(fc + l|. (2d) 

for k = 0,1,..., JV - 2 where J w = J]^ 1 is the 

identity. For N = 2 these reduce to the standard Pauli 
operators z 2 ,o — a z , j/2,0 = fj/ and a; 2 ,o = cr^. 

Lemma: Let < to < t\ < 00 and t E (to,ti) C R + . 
Both F(i) and /(t) are scalar time functions defined on 

R. If dF J? = fit) for f e (to)^i) and is a constant 
Hamiltonian, then the solution of the operator differen- 
tial equation 



we shall consider two kind of performance indices, the 



iX(t) = f(t)HX(t) 



(3) 



is given byX(t) = e -i[F(t)-F(to)]H x ^ with the initial 
state X(t ). 

Setting AF(t) = F(t) — F(to) evaluation of the ma- 
trix exponential for H — ZN,k and H = yN,k yields the 
explicit formulas 

e -iAF(t)z N , k =e -iAF(t) {7jv + [e i2AF(t) _ + + 1|} 

(4a) 

e -iAF{t)y N , k = j N k CQS _ iyN k gin 

+ I N ~lN,k- (4b) 

This Lemma implies that one have degree of freedom 
to construct control Hamiltonian to steer the quantum 
system to the target state. In this paper we focus on 
controls given by both a local sine waveform 



fu(t;t ,t u A) 



Asin^£M te[t , tl ] 
otherwise 



and a local n-order-polynomial- function fi n (t; to, ti, A) 
fin(t; t ,h,A) - 



_ A [ 2t-(ti+t ) ]n + ^ t€ [to±tl )tl ) 

otherwise 

(6) 

We further have /i s (f; i ,ii, A)d* = 

MC|z*2l ) fr Q f ln (t-t Q M,A)dt = ^3^, 
/ 00 |/i.(*;io,*i > A)| 2 df A2(t r* o) and 

jo" ^(titti.tx, A)i»«ft = ^vgr 2 - 

When some waveform controls are chosen to achieve 
a state transition at the target time tf, we still have 
some degree of freedom to optimize the control magni- 
tudes with regard to a chosen performance index. Here 



transition time Jt = ldt 
energy performance index 



tf and a combined time- 



J te = t f + X- 1 ■ / E(u(t))dt = / [1 + X~ 1 E(u(t))]dt 
Jo Jo 

(7) 

that takes into account the energy cost E(u(t)) of the 
control vector u(t) as well as the time required. For ex- 
ample, if H(t) = Yli u i{t)Hi 1 then u(t) = (ui(t)) and 
E(u(t)) = ^Z i \ui{t)\ 2 . Here A > is introduced as a ratio 
parameter that defines the relative weight of the energy 
and time resource costs, and the equivalent physical unit 



oiXisW = J ■ 



N ■ m ■ s 



(kg) ■ m 2 ■ s 3 . It 



should be emphasized that J te is reduced to Jt if A — > 00, 
so the time-energy performance Jte can be interrupted as 
the generalization of time performance. 



III. TRAJECTORY-CONSTRAINED OPTIMAL 
LOCAL SINE- WAVEFORM TRANSITION 
CONTROLS 

Consider an N— level quantum system which is gov- 
erned by the Schrodinger equation (in units of h = 1) 



N-2 

E 



= Y] [Uy,h{t)VN,k + U Z:k {t)z NM M{t)) (8) 



where yN,k and zjq,k are defined in Eq. ([2]). Such a drift- 
free Hamiltonian can be obtained for many systems by 
transforming to a rotating frame, suitably expanding the 
control fields and making certain simplifying assumptions 
as discussed e.g. in Ref. [2l| . 

It has been demonstrated in Ref. [2l| that one can 
steer the system (jHJ from an arbitrary initial state |^>o) 
to a target state \ip s ) by using (at most) AN — 5 Bang- 
Bang control based on the paramctrization of the ini- 
tial and target states in terms of the 2(N — 1) geomet- 
ric parameters (6 ; $0) = (<??, ■ • • , 0° N -x,<t>i, • • ■ , 4>°n-i) T 
and (O s ;$ s ) = (0f , . . . , 9 S N _ 1 ; </>f , . . . , (j)^^ described 
in Section 2. One can also achieve the same state tran- 
sition by using (at most) AN — 5 local sine waveform 
controls and the design process can be divided into three 
stages: (1) Steer the quNit from (0 O ;$ O ) to (0q;O) by 
performing (N — 1) local Z- rotations; (2) Transfer the 
quNit from (0 O ; 0) to (0s; 0) by performing 2N — 3 local 
Y- rotations; (3) Manipulate the quNit from (0 S ;O) to 
(0 S ; $ s ) by performing (N — 1) local Z-rotations. 

Let m = 2,3, N — 1 and j — 1, 2, N — 1, we have 

Uyfi(t) = Sgn(0° - 6 S 1 )fls(t;t 2 N-3,t2N-2,A g as) 
u y,m-l{t) — fls (t] t2N-2-m i ^JV-l-m, AqoJD N >m -l 
^ fls{t; t2N-i+rn, t2N-3+m, Affs^ )yN,m- 1 

Uz,j-l(t) = sgn(?r - <p° J )f ls (t;t J _ 1 ,t : j,A 4> a) 

= sgn(0j - 7r)/ ;s (t;t 3 Ar_5 +i ,t 3A r_4 +i , A r ) 

(9) 



3 



with control magnitudes Ago , Ags , Ago* , A+o , and Ad 

m m 1 Vj 

corresponding to the geometric parameter 0^, 0^, 0® - 



$ and 



min{0°,27T-0°}7r 
TXT I2JV-1-T 



2AT-2-TO 



2W-2 



- £ 



2N-3 



_ |g£— gfj-TT 



t 



and i 



, t2N-3+j ~ 
min{2ir — 0°,<ji°}7i 



3JV-5+J - iA- 

= tiN-5 — to and the cor- 



2JV-4+J - 

Both the transition time tf 
responding energy cost can be expressed in terms of the 
initial and target state parameters as well as the mag- 
nitudes of sine waveform functions. This implies that 
one can optimize the magnitudes of sine wave function in 
terms of time-energy performance index given by Eq. (J7J . 

Suppose the control amplitudes are bounded by L, and 
denote L* = min(L,V2X) and wi s (x) = + the 
optimal amplitudes for the sine waveform controls, min- 



imizing the performance index J te , are A^ 
Ago — Ags — Agos — L* and we have 



A 



Jt e = t f + \- 1 



E{u{t))dt = {C 1 +C 2 )w ls {L^) (10) 



where C\ and C 2 are given by 



N-x 



ci = £(*?+*f) + i*i-*;i 

2=2 

C 2 =Y^ [min(2vr -^14)+ min(27r 



(11) 



k=l 



The "best" time performance for bounded controls is ob- 
tained by letting A — > +oo. If L —> +oo, then time 
performance index J t — > 0. 

To intuitively understand the aforementioned analysis, 
we further consider a concrete example by setting N = 2, 

IV^o) = |0), \i/>.) = a^|0) + A = 2 and L = 1. For 

this concrete example, J 4 * e = and the corresponding 
constrained optimal bounded sine-waveform controls are 
given as follow 



St — 7T 



otherwise 



(12) 







otherwise 



This implies that the desired state transition is achieved 
by 2-rotations as shown in Figflja). One can also achieve 
the same state transition along one-rotation trajectory 
(FigJTJb)) by choosing 



;(*) = <(t) = 



I S i n ±VM te[0 



V2ir 2 







otht 



(13) 



and the corresponding time-energy performance is J t * e = 
3v/ ^ 7r Therefore, it is exemplified that the obtained op- 
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timal sine- waveform controls proposed in this section are 



TABLE I: Trajectory-constrained minimum transition times 
Jt for bounded controls with amplitude bounds given by L 
and optimal time-energy performance Jte for both bounded 
and unbounded controls for different waveforms. 



Case J t Jte (bounded) 

BB \U wi b (L* B ){d + C a ) 

LN !&t. w ln {L* n ){C, + C 2 ) VC"+W0A (Cl + c 2 ) 



Jte (unbounded) 
VA(Ci + C 3 ) 



LS 



[U wULDid + C* 2 ) 



2n 



(Ci + C 2 ) 



not globally optimal because they are subject to con- 
straints on the control waveforms and trajectories. It 
should be also underlined in the aforementioned example 
that one can generate the target state \ip a ) = ^|0) + 
at t k = + k = 0,1,2,...) by perform- 

ing the sine waveform functions u*(t) = u*(t) = sin 
on the controlled qubit. 




(a)2-rotation trajectories (b)l-rotation trajectory 

FIG. 1: Two kinds of trajectories on the Bloch sphere 

IV. DISCUSSION AND CONCLUSION 

Following similar analysis as in Section 3, one can 
obtain trajectory-constrained optimal nst-order polyno- 
mial waveform controls. Denote Wib(x) = |(- + f) 

w m{x) = (jg + g^ ), L* B = min(L,V\)*L* n = 
min(L, V( 2 "+iK 2 "+ 2 ) ^) and u = ( Cl + (j 2 y L with con _ 
stants C\ and C 2 defined as Eq. (fTTj) . The time and 
time-energy performances are further summered in Table 
Ufor constrained optimal Bang-Bang (BB) controls, con- 
strained optimal local sine-waveform (LS) controls, and 
constrained optimal local n-order-polynomial-function- 
waveforms(LN) controls. 

In this technical communique, we utilize the geomet- 
ric paramctrization of quantum states and the proper- 
ties of generalized Pauli operators to develop analytic 
control schemes and then construct optimal local time- 
continuous function controls to achieve state transitions 
for multi-level quantum systems. We use the remain- 



4 



ing degrees of freedom to optimize the magnitude pa- 
rameters of the time-continuous function with regard to 
time-energy performance. Constrained optimal controls 
of both sine- waveform and n-order-polynomial-function- 
waveform are obtained. When n — > oo, constrained op- 
timal LN controls approach to the constrained optimal 
BB controls, generalizing the results in (2lj . 

The choice of a time-energy performance index given 
by Eq. Q is motivated by experimental feasibility consid- 
erations, which require that the desire for fast state tran- 
sition be balanced against the need to limit the amount 
of energy required to achieve the transition. It should be 
underlined that Eq. is reduced to J t when A — > oo, 
therefore time-energy performance can be interpreted as 
a generalization of time performance. 



Although the trajectory-constrained optimal time- 
continuous waveform controls are not globally optimal, 
the resulting constrained optimal controls have the ad- 
vantage of being of a simple form, which can be given 
analytically without the burden for numerical optimiza- 
tion. It is also exemplified in section 3 that periodic 
sine waveforms controls could possibly be utilized to pe- 
riodically generate the target state, and this observation 
implies that bounded control with finite frequencies can 
be used for quantum state engineering. 
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